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ABSTRACT 

Spectropolarimetric measurements at moderate spectral resolutions are effective tracers of stellar magnetic fields and circumstellar 
environments when signal to noise ratios (SNRs) above 2000 can be achieved. The LRISp spectropolarimeter is capable of achieving 
these SNRs on faint targets with the 10m aperture of the Keck telescope, provided several instrumental artifacts can be suppressed. 
We describe here several methods to overcome instrumental error sources that are required to achieve these high SNRs on LRISp. 
We explore high SNR techniques such as defocusing and slit-stepping during integration with high spectral and spatial oversampling. 
We hnd that the instrument flexure and interference fringes introduced by the achromatic retarders create artificial signals at 0.5% 
levels in the red channel which mimic real stellar signals and limit the sensitivity and calibration stability of LRISp. Careful spectral 
extraction and data filtering algorithms can remove these error sources. For faint targets and long exposures, cosmic ray hits are 
frequent and present a major limitation to the upgraded deep depletion red-channel CCD. These must be corrected to the same high 
SNR levels, requiring careful spectral extraction using iterative filtering algorithms. We demonstrate here characterization of these 
sources of instrumental polarization artifacts and present several methods used to successfully overcome these limitations. We have 
measured the linear to circular cross-talk and find it to be roughly 5%, consistent with the known instrument limitations. We show 
spectropolarimetric signals on brown dwarfs are clearly detectable at 0.2% amplitudes with sensitivities better than 0.05% at full 
spectral sampling in atomic and molecular bands. Future LRISp users can perform high sensitivity observations with high quality 
calibration when following the described algorithms. 

Key words. Instrumentation: polarimeters - Instrumentation: detectors - Techniques: polarimetric - Techniques - spectroscopic - 
Methods: observational 


1. Introduction 


The Keck telescope on Mauna Kea and the Low-Resolution 
Imaging Spectrograph LRIS delivers a very large collecting area 

i lOm) to a moderate resoDtion Cassegrain slit spectrograph 
Goodrich & Cohenll2003l : lOke et al.lll995h . The polarimetric 
mode, called LRISp, provides a dual-beam mode with the ability 
to measur e circular and linear polarized spectra u sing achromatic 
retarders (iGoodrich et al.l[l9^ lGoodrichlll991h to modulate the 
incoming polarized light. 


In the field of stellar magnetism, several recent theoretical 
and observational studies show that resolutions of only a few 
thousand are required to detect signatures of global fields and 
star spots. Magnetic fields in TiO ba nds have been detected in 
M-dwarfs. dBerdvugina et al.ll2008bl) . Iron hydride (FeH) and 
chromium hydride (CrH) bands are also observable and mod¬ 
eled to be detectable with high polarimetric sensitivitv (Afram 
et al. I 2 OO 8 LI 2 OO 7 I : iKuzmvchov & Berdvuginall2Q13h . The 3- 
dimensional structure of star spots can be constrained w ith ob¬ 
servations in multiple molecular bands (lBerdvuginal201 ih . With 
LRISp, several scientific investigations are possible provided 
that signal-to-noise ratios (SNRs) of over 1000 can be de livered 
as we have been pursuing (IKuzmvchov et al.ll2014ll2()T^ . 


Many galactic sources are highly polarized showing de- 
tectable continuum and line pol arization effects of up to 20% 
(iTran et al.l 1199.51 : iTranIll995bHah . LRISp has been used to 
search these targets ach ieving roughly 1% polarimetric sensi- 
tivites (iTran et al.ll201 Ih . 

To achieve these high SNRs, not only are many photons re¬ 
quired, but a thorough calibration and correction of many in¬ 
strumental artifacts must be performed. Most instruments suffer 
from several problems that create artifacts in polarimetric data. 
For example, instrument flexure causes wavelength drifts of a 
fraction of a pixel between exposures or within the two polar¬ 
ized beams of a single exposure. These small instrumental wave¬ 
length shifts mimic signatures from stellar magnetic fields which 
are also w avelength shifts of sp ectral lines through the Zeeman 
effect fc.f. lBagnulo et al.l(l2013h ). A typical polarization calcula¬ 
tion requires combining spectra from several exposures. Rotat¬ 
ing retarders, as in LRISp an d other instruments, can in troduce 
several polarimetric artifacts (iHarrington & Kuhnll2008h . Instru¬ 
mental wavelength instabilities can become a serious limitation 
even with shifts as small as 0.1 pixel. The flexure of LRIS was 
measured in__2011_b^^elesco£e_staffJo_^e_ove£_a_^bce^Keck 
Telescopel201 Ih similar t o previous LRIS flexure measurements 
(ICohen & Shonhelllll996h . 
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Polarimetric instruments typically use retarders to modulate 
the incoming stellar light and a polarizer to function as an an¬ 
alyzer. Dual beam instruments utilize polarizing beamsplitters 
so that orthogonally analyzed polarization states are recorded on 
the CCD with high instrument throughput and the ability to re¬ 
move systematic errors. By differencing and / orratioing spectra 
recorded with different retarder orientations through the two or¬ 
thogonally polarized dual beams, several instrument systematic 
errors can be removed. This is typically called beam swapping. 
Several examples are included in these references: fdel Toro Ini- 


estal2003LlBagnulo et al. l2009l:lTinbergenll2005HSnik & Kelleii 

1201 sT Semel et al. 1993t Gehrelslli974 ). Dual beam instrumen¬ 
tation with beam swapping modulation techniques can cancels 
out several effects to first order. Depending on how the polar¬ 
ization spectra are computed, there are still second-order arti¬ 
facts remaining from instabilities in time (seeing, pointing jit¬ 
ter, sky transparency) and between the beams (transmissions, 
flat-helding, differential aberrations, CCD imperfections, etc). 
Since polarization spectra are computed as differences between 
measured intensities , instrumenta l artifacts must be carefully re¬ 
moved or mitigated (lKelleiill996ll . 


Several kinds of retarders often have small but detectable 
interference effects, called fringes i n some texts and are mor e 
generally polarized spectral fringes (IClarkell2005l : ISemelll2003h . 
Many retarders are manufactured as multiple layers of birefrin- 
gent materials that produce Fabry-Perot type etaloning which in¬ 
troduces fringes. These fringes are different for both the fast and 
slow axis orientations, creating spurious instrumental polariza¬ 
tion. Some of the super-achromatic type retarders with multiple 
layers can have fringes producing spectral intensity modulation 
of well over 1% amplitudes. These fringes are often modeled or 
removed to some residual err or level with various fun ction fits or 
Fourier filtering techniques. dAitken & Houg h||200l|; Handes & 
Howarth i 19961 Adamson & Whitted 1 1 9951 : IClarkel120051 : ISemell 
I2003h . 


Telescope and instrumental polarization is important because 
knowledge of the continuum polarization can be used in addi¬ 
tion to line polarization as a constraint on circumstellar environ¬ 
ments. Because this dual-beam spectropolarimeter is mounted 
at Cassegrain focus in a mostly symmetric optical beam, it 
has quite minimal instrumental polarization. The instrument is 
mounted roughly 10 arc minutes off-axis showing low polariza¬ 
tion induced by the telescope. Most telescopes even in axially 
symmetric beams show instrumental induced polarization at the 
0 .1% level from asymmetries in the optical coatings, oxidation, 
metallic properties, etc. These small but significant telescope po¬ 
larizations are seen in several imaging instru ments such as Plan- 


iHough et al.l 20061; Bailev et al 

l2008l 20101; Wiktorowicz & 

Matthews l2008t iBerdvugin et al. 

20061; iBerdvugina et al.ll201 iL 


data reduction artifacts can all introduce continuum polarization. 


Spectrographs generally are limited in their stability by 
changes in the optical path (pointing, slit tracking jitter, flexure, 
dispersive optics sensitivities). Spurious instrumental polariza¬ 
tion can also be caused by incomplete scattered light compen¬ 
sation, CCD instabilities, polarization induced by reflective op¬ 
tics (e.g. oblique fold mitTors) and imperfect coatings on optics. 
Thus, measuring the absolute value of the polarization at high 
accuracy across the entire continuum of a dispersed spectrum 
presents challenges to both the instrument and the data analy¬ 
sis pipeline. Many spectrographs have much higher polarization 
sensitivity across individual spectral lines because the contin¬ 


uum variations can b e differentially subtr acted across a small 
wavelength region (cf. lPerevra et alJ (l2015h l 

In addition to polarization being created by the telescope, 
the instrument can also scramble or mix incoming polarization 
states. This mixing, called cross-talk can be from the unpolar¬ 
ized intensity to detected linear or circular polarization or simply 
mixing between linear and circular states. The LRISp retarders 
are highly achromati c, reducing th e mixing between linear and 
circular polarization (ICohenll2005h . 

The LRISp instrument was upgraded to include a second 
blue camera which also,can be used in nolarimetric mode (Co¬ 
hen '2005|;lMcCardiiLeLalJLl998D. In 2009, the LRIS red detec¬ 
tor was upgraded to include higher sensitivity at longer wave¬ 
lengths (iRockosi et alJl2010r) . The atmospheric dispersion cor¬ 
rector ( ADC) was mounted in 2007 and includes transmissive 
prisms dPhillins et al.ll206^ . This ADC was tested in 2007 to 
only marginally impact the measured d egree of polariza tion and 
angle of polarization for standard stars (lGoodrichll2007t) . 


2. Observations 


We observed a range of targets on August 22nd and 23rd 2012. 
The 831/8200 grating was used for the red channel at an angle of 
37.47° giving coverage from 789nm to 1026nm. The blue chan¬ 
nel used the 300/5000 grism with the 680 dichroic with reason¬ 
able sensitivity from 380nm to 776nm though drastic throughput 
losses were seen long of the 680nm dichroic cutoff. 

The target list included magnetic stars, brown dwarfs and 
a range of calibration standards. EV Lac and V1054 Oph are 
magnetic flare stars of roughly M3 to M4 type. These stars have 
roughly known magnetic field stren gths and have been studied 
exten s ively in the optical and radio ( Pettersen et al.l[T984l : ISaaij 
1 19941 : Llohns-Krull & Valentil 1 19961 l2000ll . We use them here 
as stars where we expect relatively large and detectable signals. 
The star HD20630 is a G5Vv star of BY Dra type. Though this 
star is magnetic with known variability, it is listed as a bright 
unpolarized standard star (in continuum filters) on the UKIRT 
standard star listQ and in several publications (lGehrelsl|1974ll . 
Ceres is the largest main belt asteroid and has a visual magnitude 
ofV=8. 


2.1. Polarimetric Modulation and Demodulation 


In this paper we use the standard Stokes vector formalism to 
describe polarized spectra. Linear polarization is denoted as Q 
and U while circular polarization is V. When we normalize a 
spectrum by the total intensity, we use lower case symbols. For 
instance q - Qjl. 

We use the general framework for measuring polarization as 
a modulation and demodulation process. The Stokes parameters 
are typically described as differences between intensities mea¬ 
sured with retarders at different orientations. In the limit of per¬ 
fect instrumentation and achromatic optics, a spectropolarimeter 
can create exposures that mimic the definition of the Stokes pa¬ 
rameters. Several modulation strategies are in use in solar, space 
and night time applications in order to balance the need for ef¬ 
ficiency, redundancy, error ch ecking through null spectra and 


2003; 

del Toro Iniesta & ColladosI 20001; ISnik & Kellei 

2013; 

Nagaraiu et alJ 12007 

; iTomczvk et al.l 120101: ISnik et al, 

I2OI2I 


2009I : Ide Wiin et al 


201 oh . In the typical notation, the i 


instru¬ 


ment modulates the incoming polarization information in to a 


http://www.jach.hawaii.edu/UKIRT/instruments/irpol/irpol_stds.html 
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Table 1. Observed targets for LRISp August 22nd and 23rd. 


Name 

Date 

Exp 

Viv 

Elevation 

Defocus? 

Spec 

Type 

EVLac 

22 nd 

20 

2000 

64-65 

No 

M4.5V 

flare star 

V 1054 Oph 

22 nd 

15 

2000 

59-60 

No 

M3.5Ve 

flare star 

HD 20630 

22 nd 

1.2 

1000 

73-74 

Yes-H 

G5Vv 

star type BY Dra 

2MASS 

22 nd 

600 

350 

67-84 

No 

L3.5 

Brown Dwarf 

2MASS 

22 nd 

600 

350 

88-77 

No 

L3.5 

Brown Dwarf 

2MASS 

22 nd 

600 

350 

72-59 

No 

L3.5 

Brown Dwarf 

LSRJ 

22 nd 

600 

1000 

74-77 

No 

M8.5V 

Brown Dwarf 

LSRJ 

22 nd 

600 

1000 

74-64 

No 

M8.5V 

Brown Dwarf 

LSRJ 

22 nd 

600 

1000 

62-49 

No 

M8.5V 

Brown Dwarf 

Ceres 

23rd 

60 

2000 

59-54 

Yes- 

Asteroid 

Solar 

HD 20630 

22 nd 

2 

2000 

73-74 

Yes 

G5Vv 

star type BY Dra 

HD 174160 

22 nd 

5 

1200 

73-70 

Yes-H 

F8V 

star 

2MASS 

22 nd 

600 

200 

63-77 

No 

L3.5 

Brown Dwarf 

2MASS 

22 nd 

600 

200 

81-88 

No 

L3.5 

Brown Dwarf 

2MASS 

22 nd 

600 

200 

79-64 

No 

L3.5 

Brown Dwarf 

LSRJ 

22 nd 

600 

1200 

77-74 

No 

M8.5V 

Brown Dwarf 

LSRJ 

22 nd 

600 

1200 

73-63 

No 

M8.5V 

Brown Dwarf 

LSRJ 

22 nd 

600 

1200 

59-46 

No 

M8.5V 

Brown Dwarf 


This table shows all complete polarimetric data sets used for development of this data reduction pipeline. Note that the star we de¬ 
note as LSRJ is LSR J18353790+3259545 and the star we denote as 2MASS is an L3.5 brown dwarf: 2MASS J00361617+1821104. 
The star HD 20630 is a magnetic star (BY Dra type) but it is als o listed on the UKIRT IRPOL list of unpolarized standards 
http://www.jach.hawaii.edu/UKIRT/instruments/irpol/irpol_stds.html (lGehrelslll974h . The statistical upper limit to the signal-to- 
noise ratio of each polarimetric exposure was computed empirically from polarimetric data. We computed the pixel-to-pixel 
variance in the polarization spectra after applying high-pass filters to isolate the statistical noise. This noise limit represents the Viv 
limitations from the detected photon flux as an upper limit to the data sensitivity. We also show the azimuth and elevation range for 
the telescope to illustrate the variation in local gravity during an exposure. The slit de-rotator was used and set to parallactic for 
each brown dwarf exposure adding to the gravitational orientation changes between exposures. The telescope was defocused for 
several bright targets to investigate this technique for increasing the exposure time to saturation. The Yes- indicates some defocus 
while YesH- indicates substantial defocus for very bright targets. The spectral classification and star type from SIMBAD are shown 
in the last two columns. See the text for details. 


series of measured intensities (I,) for i independent observations 
via the modulation matrix (O,;) for j input Stokes parameters 

(S;): 


1 / = O^jSj 


( 1 ) 


In most night-time polarimeters, instruments and associated 
data analysis packages a modulation matrix that separates and 
measures individual parameters of the Stokes vector as well as 
providing redundan t information for use in characterizing ins tru- 
ment performance (ISemel et alJll993L iDonati et alJll997h . In 
the Stokes definition modulation scheme, there are 6 exposures 
recorded each corresponding to an independent Stokes parame¬ 
ter (QUV). 


( 


0,v = 


V 


1 

1 

1 

1 

1 

1 


-Hi 0 Of 

-10 0 

0 -Hi 0 

0-1 0 
0 0 -Hi 

0 0 -1 j 


( 2 ) 


In ESPaDOnS, FORS and other instruments, additional re¬ 
dundancy is achieved by making another set of measurements 
using the same mod ulation matrix but w i th all retarders ro- 
tated by 180 degrees dBagnulo et alJl2()()^ iDonati et al.lll999l : 


ISemel et'anil993h . In LRISp, this type of modulation is accom¬ 
plished in two separate optical configurations (for two separate 
exposures). A rotating half-wave super achromatic retarder plate 
(HWP) is mounted in front of the analyzer. The HWP is ro¬ 
tated in a sequence of [0°, 45°, 22.5°, 67.5°] in order to accom¬ 
plish linear polarization modulation in 4 exposures with beam 
swapping. Circular polarization is measured by rotating a sec¬ 
ond quarter-wave achromatic retarder plate (QWP) into the opti¬ 
cal path using the calibration filter wheel. This QWP is fixed in 
a single orientation and modulation is accomplished by rotating 
the HWP by 0° and 45° behind the QWP. 

In the Stokes definition scheme, calculation of each Stokes 
parameter from intensity spectra follows Equation [3] is imple¬ 
mented as a series of normalized intensity differences recorded 
in two exposures assuming perfect modulation and achromatic 
optics: 


<2 

— ^qo + q\ 


h - h _ h - h 
h) + h h + h 


(3) 


We wish to highlight that these normalized intensity differ¬ 
ences when assumed to represent a Stokes parameter can in¬ 
troduce several types of instrumental errors while also ignor¬ 
ing cross-talk. In dual beam systems, two pairs of spectra are 
recorded in two exposures. Thus each part of the ratio is subject 
to instrumental uncertainties that are introduced between expo¬ 
sures. 
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In many night time spectropolarimeters, the instrument is de¬ 
signed so the cross-talk is below some nominal design value. 
Chromatic effects are often just minimized and subsequently left 
uncalibrated, but no additional spectra are used to compute a 
Stokes parameter. By using additional spectra, or additional cal¬ 
ibrations, cross-talk can be further minimized. Each source of 
polarimetric error must be considered when choosing an optimal 
modulation scheme and associated data processing algorithms. 

Other instruments choose to modulate and measure all 
Stokes parameters using only four mea surements or pursue less 
redundant but more effic^nt»chemesjldeJ^ijne^a^_2Ql(]Ldel 


2012 

2009 

ilNaearaiu et al.ll2007; Keil et al. 2011; Elmore et al. 

201C 

1992 

). Eor instance, one can use alternate retardance and 


fast axis orientations to give a modulation matrix that uses only 
four exposures to measure a Stokes vector with maximal effi¬ 
ciency. 

One recovers the input Stokes vector from the series of in¬ 
tensity measurements by inverting the modulation matrix (O). If 
the matrix is not square, and non-unitary (as in the Stokes def¬ 
inition scheme) one can simply solve the over-specified system 
of equations via the normal least squares formalism: 


S = 


0^1 

0^0 


(4) 


Other modulation schemes are easily crated using tunable 
liquid crystals and the modulation matrix does not need to have 
any particular symmetry. We have implemented this ours elves on 
other spectropolarimeters dHarrington et al.ll2010l.l20llT) . In sys¬ 
tems with more complex or less redundant modulation schemes, 
additional calibrations with upstream polarizers and retarders are 
often used to achieve the highest calibration accuracy and re¬ 
move residual chromatic effects from the demodulation process. 
The demodulation process of Equation |4] can be used regardless 
of modulation scheme. 

Though often not performed, these same kind of demodu¬ 
lation processes can be applied to Stokes definition type mod¬ 
ulation schemes to remove residual chromatic errors if the im¬ 
perfections are estimated through a calibration procedure. Eor 
instance, in LRISp, there are calibration polarizers mounted in 
the filter wheel ahead of the rotating HWP retarder. This polar¬ 
izer can be used to measure some of the chromatic properties of 
the HWP and to modify the modulation matrix of Equationl^to 
account for imperfections. Polarized standard stars or more elab¬ 
orate calibration optics can be used to derive the system Mueller 
matrix to correct for some residual cross-talk. Alternatively, the 
daytime sky calibrations we outline here and elsewhere can be 
used to measure the system Mueller matrix and to apply correc¬ 
tions to the demodulated spectra to acc ount for any uncalibrate d 
cross-talk to some residual error levels dHarrington et al.ir201 ih . 
Harrington et al. 2015. 


2.2. Geometric calibrations - wavelength stabilization 

The hrst steps in spectral extraction is locating the spectral orders 
and identifying the basic optical configuration. When combining 
6 exposures to make a single quv spectral data set, any instru¬ 
mental instability can produce spurious signals. In general, spec¬ 
tral order curvature, anamorphic magnihcation, and tilt of the slit 
image against the ccd pixel grid all can impact polarimetric data. 
This is especially true given telescope guiding imperfections and 
instrument flexure. Eor the LRISp red channel, the basic param¬ 
eters of our optical extraction are shown in Eigure [T] We derive 


several parameters from fits to several arc lamp calibration expo¬ 
sures. The spectral resolution is derived from arc lamp calibra¬ 
tion exposures. The resolving power (R=/i / SA) shown in Eigure 
[T]is R=2500 at SOOnm rising to R=3500 at lOOOnm. In this con- 
hguration, the spectra are oversampled. Erom Gaussian fits to arc 
lamp spectral lines, we hnd the spectra to be sampled at 4.5 pix¬ 
els to 5.5 pixels in the Gaussian full-width half-max (EWHM). 
With this over sampling of 0.56A to 0.59A per pixel, we can test 
for several instrumental artifacts and apply several types of data 
post processing biters to remove noise sources. 

With the arc line exposures, we measure how the wavelength 
coordinates change as the HWP rotates through the typical mod¬ 
ulation sequence of 0°, 45°, 22.5°, 67.5°. This rotating HWP 
causes a drift of roughly 0.15 pixels between the two separate 
modulation states. Eor the LRISp HWP as mounted, the offsets 
average [0, 0.08, 0.13, 0.10] pixels referenced to the hrst ex¬ 
posure. There is also a mild wavelength dependence across the 
CCD. Note that the two Stokes q exposures would show a wave¬ 
length drift of 0.1 pixels in between modulated images result¬ 
ing in imperfect subtraction introducing artifacts resembling the 
derivative of the intensity proble with wavelength. However, the 
Stokes u exposures would show substantially less wavelength 
drift between modulation states, but would be offset from Stokes 
q spectra by 0.1 pixels in wavelength. 


Spectral Pixel 

0 1000 2000 3000 4000 



800 850 900 950 1000 

Wavelength (nm) 


Fig. 1. The derived spectral sampling and spectral resolution for the 
LRISp red channel. The black parabolic curve shows the spectral sam¬ 
pling (pico meters per pixel) using the left hand y axis. The mapping 
between the wavelength solution and each spectral pixel is found by 
comparing the top and bottom x axes. The spectral resolution is derived 
as the FWHM of the arc lamp Gaussian fits as defined in the text. The 
arc lines typically have a FWHM of 4.5 to 5.5 spectral pixels. The spec¬ 
tral resolution is derived as the wavelength decided by the FWHM of 
the arc line fits and is shown with the symbols using the right hand y 
axis. The resolution was between 2500 and 3500 across the sampled 
wavelengths. The blue symbols show the spectral resolution of the top 
polarized beam. The red symbols show the spectral resolution of bot¬ 
tom polarizers beam. The red curve in between these symbols shows 
the average spectral resolution of both top and bottom beams. 


Slit guiding for Keck is software-referenced and the user can 
vary the stellar location along the length of the slit. In addition, 
we describe later how we had substantial guiding drifts tracking 
our targets. We observed at higher elevations and see expected 
drifts for an altitude-azimuth telescope with a low-bandwidth 
guider control system. With drift of the optical beam along the 
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slit, uncorrected geometrical tilt of the dispersed spectra will 
lead to wavelength drifts between exposures. 

For the red channel, we find significant tilt of the monochro¬ 
matic slit images against the pixel rows. There is roughly 2 pix¬ 
els of wavelength change on the CCD pixels from the bottom to 
the top of the imaged slit over the 300 spatial pixels sampled. 
We observed up to 40 spatial pixels of guiding drift during our 
observing run. This spatial drift combined with the spectral tilt 
would give wavelength instabilities of up to half a pixel if this 
geometrical effect is not compensated in the pipeline. 



0 1000 2000 3000 4000 

Spectrol Pixel 


Fig. 2. An example of an extracted 2-D spectrum for a typical LSRJ 
brown dwarf observation. The two beams are shown after tilt-correction 
as observed August 22nd after geometric extraction. Several night sky 
glow lines are visible as the vertical stripes. The top beam has notice¬ 
ably higher throughput than the bottom beam as was typical for this 
spectrograph configuration. The intensity was linearly scaled with blue 
and black colors corresponding to the highest intensities. 


To compensate for this geometrical tilt, the data is linearly 
up-sampled to a 0.01 pixel grid and then shifted spectrally to 
compensate for the tilt. The data is then averaged back down 
to nominal 1-pixel sampling. An example of an extraction after 
tilt correction is shown in Figure |2] In this Figure, the stellar 
spectrum is in the center of the two extracted beams. The 4000 
spectral pixels of the detector are shown on the X-axis. We only 
show a subset of the extracted spatial pixels to clearly show the 
stellar spectrum. 


3. Automatically Removing Cosmic Rays With 
Iterative Filters 

The LRIS red channel pixels are 300 microns deep, giving a 
fairl y high rate of cosmic rav hits in long exposures ( Rockosi 
et al. l2010h . An example of the background region in a typical 
10 minute brown dwarf exposure is shown in Figure |3] Bright 
night sky glow lines are seen in addition to many cosmic ray 
hits. Cosmic ray damage in this CCD often spans 10s of pixels 
and most spectral pixels are contaminated at some level, making 
sensitive polarimetry difficult. 

We define the spatial profile as a trace through the data that 
is orthogonal to the wavelength direction after geometric cali¬ 
bration of the extracted spectral data has been performed. This 
spatial profile contains the atmospheric seeing, telescope jitter, 
guiding imperfections, optical imperfections (ghosts) and can be 
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Spectral Pixel 



Fig. 3. This Figure shows a small region of the 2-D spectrum used to 
extract and remove the sky-glow lines in a single polarimetric exposure. 
The two polarized beams are shown separated by a white line. Typical 
cosmic ray hit rates are seen. 



Spotiol Pixel 


Fig. 4. An example of the cosmic ray iterative filtering process. The 
black curve shows the raw stellar spatial profile (a cut through the data 
orthogonal to the spectral direction representing the local seeing, tele¬ 
scope jitter and optical imperfections). A very large cosmic ray hit is 
seen damaging several pixels in the middle of the spatial profile at pix¬ 
els 5 to 10. The wide spread damage requires the iterative solution to 
reject the damaged data. After iterating, 5 spatial pixels are rejected and 
ignored in the shift-n-scale profile fit. The blue curve shows the median 
spatial profile used in the fit. The red curve shows the spatial profile 
replacements. The right hand y axis shows the noise estimates used in 
the filter. The triangle symbols below show the noise estimates on a per 
pixel basis used in the rejection. A 2-sigma filter was applied in this 
case. 


used to apply data-derived filters for various error sources. An 
example spatial profile is shown in Figure |4l 

An iterative method has been developed to use the spatial 
profile to filter out these cosmic ray hits based on what’s called 
optimal extraction (iHornell 19861 : iMars hl ll989ll . The spatial pro¬ 
file is computed over a range of wavelengths. We find that 100 
spectral pixels is a good compromise between increasing the 
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SNR of the spatial profile and ensuring that the wavelength of 
the spatial prohle matches that of the wavelength to be filtered. 

The first step in the hlter is to shift and scale the median 
spatial profile to the individual spatial prohle of interest. The 
least squares solution implemented over the i spatial pixels that 
independently contribute to the problem. 

We model the shift-and-scale problem with a general nota¬ 
tion where D denotes the data to be ht and P denotes the prohle 
used to do the htting. The derivative in the spatial direction is dx 
The data (D) is modeled as a sum of three terms; a constant (a) 
times the prohle (P), a constant (b) times the derivative of the 
prohle, and an additive constant (c). 


D - aP + b— + c (5) 

ox 


count level an order of magnitude larger than the detected stel¬ 
lar Hux. The hlter iterates through the shift-n-scale process and 
identihes spatial pixels with noise limits above the user-dehned 
thresholds for read and shot noise. After all points have been 
rejected, the red curve shows the values used to repair the spa¬ 
tial prohle. The triangle symbols in Figure |4] show the residual 
noise levels for the points included in the ht. The user specihed 
threshold was set at 2cr for this example. 

Typical performance of the cosmic ray hlter is shown in Fig¬ 
ure |5] The black curve shows the quv spectra with substantial 
cosmic ray hits contaminating a large fraction of the wavelengths 
covered. The blue curve shows the corresponding hltered quv 
spectra with effective removal of the cosmic ray hits. 

4. Spatial profiles and slit guider tracking 


The total error to be minimized is a sum over all spatial pix¬ 
els: 


-aP-b—- cf (6) 

ox 

I I 

In order to solve the least squares problem, we take the par¬ 
tial derivatives with respect to the three coefhcients (a,b,c) and 
set them equal to zero: 



dE _dE _ dE 
da db dc 


(7) 


After doing the partial derivatives, collecting terms and solv¬ 
ing for 0, we get a system of three equations for three variables; 
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With the least-squares solution, we can compute the differ¬ 
ence between the shift-n-scaled median spatial profile and the 
individual spatial profile to be corrected. The difference between 
these two spatial profiles is then compared against a measure of 
the expected noise at every spatial location. The expected noise 
has two contributing terms. One term is a constant background 
noise estimate representing the read, dark and scattered light 
background variation. The second term accounts for the spatially 
varying flux levels and is proportional to the square root of the 
number of detected counts. This shot-noise term is computed us¬ 
ing the shift-n-scaled median spatial profile. 

This process is iterated until convergence is achieved and no 
spatial pixels show noise levels above a user-determined set of 
thresholds. At each step of the loop, the single worst offending 
point above the noise threshold is rejected from consideration by 
setting the weighting in the ht to zero. Once the iteration has con¬ 
verged, it delivers a shift-n-scaled median prohle that hts all non- 
rejected points below the noise threshold. After convergence, all 
rejected spatial prohle points are replaced with the shift-n-scaled 
median spatial prohle values. This iteration ensures that very 
large cosmic ray strikes are properly corrected without undue 
inhuence in the htting process. 

An example of this process is shown in Figure |4l In this ex¬ 
ample, a spatial prohle for a 2MASS brown dwarf exposure on 
August 22nd is shown in the black curve. A very large cosmic 
ray hit is seen at spatial pixels 5 to 10. This cosmic ray shows a 


The Keck slit guider was used to acquire and track our targets. 
In some cases the guiding delivered spatial prohles which re¬ 
mained centered to within 5 spatial pixels for the duration of a 
6-exposure polarimetric data set. For other exposures, there was 
drift of over 50 spatial pixels. This corresponds to roughly 20% 
of the slit length (290 spatial pixels in our extractions). Thus 
our pipeline was conhgured to compute the center of light for 
each exposure of a polarimetric data set. The data was extracted 
around this detected center-of-light for each exposure to com¬ 
pensate for the guiding drift. 

A common technique to increase the exposure time until sat¬ 
uration and increase the duty-cycle of measurements is to de¬ 
focus the telescope. We tested the effects of defocus on achiev¬ 
ing high SNRs without saturation on very bright targets. During 
spectral extraction, the spatial prohle width varied from roughly 
10 pixels full width at 10% max for in-focus brown dwarf targets. 
For the bright unpolarized standards, this spatial width increased 
to over 40 pixels at 10% max. This focus shift changes the op¬ 
tical beam footprint through both QWP and HWP retarders, the 
polarizing beam splitter and all downstrea m spectrograph opt ics 
and care must be taken with calibration (cf. iTinbergenI (I20()% 1 . 

Figure |6] shows the spatial prohles for four stars to illustrate 
the change in telescope focus. The top panels show highly de- 
focused observations of bright stars. The bottom panels show an 
in-focus brown dwarf on the bottom left and a mildly defocused 
Ceres exposure on the bottom right. The data reduction pipeline 
allows for a variable spatial extraction width that changed be¬ 
tween 30 and 80 pixels to accommodate this wide ranging de¬ 
focus. In addition, by including the minimal number of spatial 
pixels needed to capture the delivered target hux, we minimize 
the impact of cosmic rays and other detector noise contributions 
(cosmetics, bad columns, etc). 

5. Flexure Compensation 

There is substantial wavelength drift even within a polarimetric 
data set due to instrument hexure. These wavelength instabilities 
cause a major source of spurious instrumental polarization. Mea¬ 
surement, compensation and accurate error budgeting of these 
types of systematic effects are critical to interpr etation of spec- 
tropolarimetric results (c.f. lBagnulo et al](l2013h l. Some instru¬ 
ments use a more redundant modulation scheme where there are 
4 exposures per Stokes parameter (using the Stokes dehnition 
modulation scheme). These redundant schemes do provide for 
error assessment through the so-called null-spectrum in addition 
to further removal of some instrumental error sources. However, 
requiring 4 exposures also introduces the possibility of wave¬ 
length jitter over long exposure times with gravity changes (for 
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Fig. 5. This Figure shows 2MASS set 3 from August 23rd. Stokes quv polarization was computed using the standard Stokes definition modulation 
assumption with beam swapping as The cosmic ray rejection filter was applied with the noise threshold of l.Scr. A 10-count background 

noise level was set. A 50-pixel smoothing width was used to compute the median spatial profile in the in the shift-n-scale fitting algorithm. The 
black lines show the computed polarization spectra before cosmic ray filtering. The blue line shows the resulting polarization after the iterative 
cosmic ray filtering is applied. Cosmic rays are present in a substantial fraction of each quv spectrum. The lower right panel shows the intensity 
spectrum with a black line in units of 100s of detected counts per spectrum. The sky glow spectrum extracted from spatial pixels outside the region 
illuminated by the star is shown in red. For 2MASS, the sky glow lines often are brighter than the star at certain wavelengths. 


Cassegrain instruments like LRISp). iBagnulo et alJ (1201 3h also 
found that using 4 dual-beam exposures (instead of 2 in our 
scheme) introduced more systematic errors through flexure in¬ 
stabilities. Additionally, some of our science targets have fast 
rotation periods and there is noticeable change over even the 2 
polarimetric exposures of a single Stokes parameter spectrum. 
Thus any science campaign must balance the error budgeting 
between the different types of systematic errors based on the 
specific use case (exposure time, flexure predictions, null spec¬ 
trum requirements, overly redundant modulation for error sup¬ 
pression, etc). 

In the brown dwarf data sets, a polarimetric observation can 
last over an hour. For faint targets, we can use sky glow lines to 
clearly align the wavelengths for each portion of the spectrum 
with a sufficient number of lines. However, in brighter targets 
with shorter exposure times, we must use the telluric absorption 
lines as an absolute wavelength reference. The 930nm to 940nm 
bandpass has many absorption lines that can easily be used for 
correlation. Figure [T] shows the telluric wavelength region for 
the first data set recorded on August 22nd. The variation in de¬ 
tected intensity is due mainly to telescope guiding drifts. The 
two different orthogonally polarized beams produced by the po¬ 
larizing beam splitter have noticeable throughput variations. The 


red curves from polarization state 0 are roughly half the flux of 
the polarization state 1 beams. 

We derive a wavelength correction for each exposure by run¬ 
ning a cross correlation analysis. All spectra are up-sampled by 
a factor of 100 to derive correlations in 0.01 pixel bins. The po¬ 
larization state 0 beam of the first exposure is used as a reference 
for all subsequent exposures and polarization states in the 6 ex¬ 
posure set. From these correlation functions we find the peak 
and derive a wavelength shift to align all exposures. These tel¬ 
luric corrections were typically less than a pixel but with varying 
behavior with exposures and between the dual orthogonally po¬ 
larized beams. 

Flexure that is common to both dual orthogonally polarized 
beams will be removed to first order with the beam swapping 
applied in the nominal Stokes definition modulation scheme. For 
the Ceres observations we use as a fringe standard, we find a 0.5 
pixel drift in wavelength but this shift is consistent between the 
two polarization states. Provided the wavelength drift between 
the two polarized beams is consistent, the primary error propor¬ 
tional to intensity derivatives will be canceled when calculating 
the polarization spectra using the difference method (a-b)/(a-i-b). 
Additional errors proportional to the second derivative of the in- 
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Fig. 6. The spatial profile of the extracted spectra after spatial trimming for guider correction and optimal extraction filtering. The unpolarized 
standard HD 174160 is in the upper left from August 23rd with an extraction width of 60. The unpolarized but magnetic BY Dra star HD 20630 is 
in the upper right from August 22nd with an extraction width of 80. The magnetic flare star EV Lac is in the lower left from August 22nd with an 
extraction width of 30. The asteroid Ceres is in the lower right from August 23rd with an extraction width of 50. All 12 spatial profiles are shown 
(6 exposures for quv and 2 polarized beams per exposure). 



Fig. 7. The intensity spectra for the LSRJ data set 1 from August 22nd. 
The detected flux is shown for each of the 6 exposures and 2 polar¬ 
ization states (dual beam). The polarization state 0 beam is red while 
polarization state 1 is blue. 


tensity will be come significant if the drift between exposures is 
large and uncompensated. 

Even though these wavelength shifts are small, they can have 
large impact on the computed polarization if differential effects 
occur. For some of our data sets, we measure wavelength offsets 
of roughly 0.2 pixels between the dual orthogonally polarized 
beams. As an example, the derived spectral pixel shifts for the 
second August 22nd LSRJ data set is shown in Figure^ There is 
a noticeable difference in behavior between the polarized beams 
for the final exposure 5. 

This difference of 0.2 pixels in wavelength for one exposure, 
though small, introduces a very large systematic error in polar¬ 
ization. Since the computed quv profiles are differences in in¬ 
tensities, any wavelength drift imprints a quv signature that is 
proportional to the derivative of the intensity with wavelength. 
This small 0.2 pixel wavelength drift is enough to imprint a 
0.5% Stokes V signature change when the corresponding inten¬ 
sity spectra has a strong absorption line. Applying a wavelength 
drift correction is critical to deriving accurate polarization spec¬ 
tra with LRISp. 

Another added benefit of including cross-correlation analy¬ 
sis in our LRISp pipeline is that all spectra can be easily refer¬ 
enced to a common wavelength. Figurej^shows the wavelength 
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Fig. 8. The pixel shifts for the LSRJ data set 2 from August 22nd. The 
polarization state 0 beam is red while polarization state 1 is blue. 


shift (in pixels) derived from correlating the Ceres telluric spec¬ 
tra with all exposures for all August 22nd and 23rd data sets. The 
flexure derived roughly +2 pixels within a night and a fraction 
of a pixel within a data set. 
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Fig. 9. The wavelength drift (in pixels) for all quv data sets derived 
through cross correlation of telluric lines with the nominal Ceres data 
set on August 23rd. The symbol for q is the triangle, u is the diamond 
and V is the cross. The left side shows August 22nd data while the right 
side shows August 23rd. The wavelength drift is roughly ± 2 pixels. 


Additionally, we have tested the wavelength stability across 
the detector. In principal, the flexure could induce wavelength 
changes that are not completely removed by a single shift for all 
wavelengths. Any optical distortion or second order effects could 
create a more complex functional dependence on the wavelength 
solution with flexure. We ran a cross-correlation of the intensity 
spectra within a complete polarimetric data set. Essentially we 
find that the wavelength solution is well corrected by a single 
shift of all spectral pixels to within at least 0.05 pixels. We And 
perturbations in wavelength regions where we detect spectropo- 
larimetric signatures, as we would expect from magnetic held 
effects. As expected, with the highly defocused standard stars 
without the use of the ADC, there are shifts in the wavelength 
solution of up to 1 pixel across the detector. However, this drift 
only influences the unpolarized standard star observations. 


6. Spectral Fringes 

Another feature of LRISp data at high SNR is a spectral fringes 
caused by interference within the achromatic retarders. Fringes 
such as these can be common in night-time spectropolarimeters 
such as the Intermediate dispersion Spectrograph and Imaging 
System (ISIS) on the William Herschel Telescope (WHT) or the 
Anglo-Australian Telescope (AA T) with the Royal Green wich 
Observatory ( RGOl spectrograph (lAitken & Hou gh||2001|; Har- 
ries & Howarth fl 996l : lDonati et al.l 1 999h . For most of these spec¬ 
trographs over relatively limited wavelength regions, the spec- 
tropolarimetric fringe follows a simple functional form. Chirp 
functions or simple harmonic Alters in Fourier space are enough 
to suppress the fringes below the statistical noise limits. 

For LRISp, the fringe depends on the Stokes parameter be¬ 
ing measured. Stokes qu measurements use only a single rotating 
half-wave retarder (HWP) giving a fringe at 0.2% amplitude at 
850nm. Measuring Stokes v with LRISp involves adding a sec¬ 
ond achromatic quarter wave plate (QWP) fixed in front of the 
HWP. This second optic complicates the fringe and increases the 
amplitude to roughly 0.5% at 850nm. Figure [TOl 
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Fig. 10. This Figure shows the quv spectra in the 850nm region from 
EV Lac. The SNR is estimated to be about 1200 to 1800 in each Stokes 
parameter. A strong spectral fringe is seen dominating the quv spectra. 
The fringe amplitude is well above the statistical noise limits and is well 
sampled with many spectral pixels. 


The exact form of the fringe power spectra depend on several 
things. First, the telescope was defocused for the brighter targets. 
During spectral extraction, the spatial profile width varied from 
roughly 10 pixels full width at 10% max to over 40 pixels. This 
means the beam footprint as the light passed through both QWP 
and HWP retarders was substantially different. Figure |6] showed 
the spatial profiles for four stars to illustrate the change in tele¬ 
scope focus. This defocus changes the size of the beam and also 
where the beam from each held angle passes through the optic. 

There are two common methods to subtract this fringing. 
First, an unpolarized standard star is observed under the same 
conditions and the corresponding spectral fringe is subtracted 
from the science target spectra. Second, Fourier Alters can be 
applied to the science target spectra without any need for a cal¬ 
ibration target provided the frequencies to be Altered are known 
for the specific instrument configuration. 

In the first method, the calibration standard can be observed 
with very high SNRs and simply subtracted without degrading 
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the SNR of the target. This assumes that the spectral fringe is 
stable between calibration standard and science target. One po¬ 
tential error is that telescope guiding drifts can change the angle 
and location of the beam through the optic. Another is that grav¬ 
ity or time dependent instrument changes (flexure, temperature, 
focus, etc) can substantially alter interference fringes. 



Wdvtitrtgih (mn) 

Fig. 11. The spectral fringes present in 4 observed calibrator targets 
for the 855nm to 865nm wavelength region. The unpolarized standard 
HD 174160 from August 23rd had a SNR of 1200, the unpolarized but 
magnetic HD 20630 star from August 22nd had a SNR of 1000 and 
August 23rd had an SNR of 2000. Ceres, though presenting continuum 
polarization, matches the unpolarized standards from August 23rd with 
an SNR of 2000. The blue curve shows the average of all 4 individual 
spectra. 


Figure[TT]shows the quv spectra for the four possible calibra¬ 
tion standards we observed August 22nd and 23rd. A clear and 
repeatable spectral fringe signature present. Ceres had a of SNR 
of 2000. The two unpolarized but magnetic HD 20630 observa¬ 
tions had SNRs of 1000 on August 22nd and 2000 on August 
23rd. The unpolarized standard HD 174160 had a SNR of 1200. 
Though Ceres is an asteroid with a continuum polarization, it is 
not expected to show significant spectropolarimetric signatures 
in this wavelength region. Since the telescope focused changed 
substantially between the unpolarized standards and Ceres, some 
systematic variation is expected between the observations. There 
is some small but statistically significant variation between the 
fringes detected in Figure[TT] We note that the Ceres spectra are 
very similar to all the unpolarized standard observations. Addi¬ 
tionally, the Ceres spectra were only mildly defocused and match 
the brown dwarf focus position much closer than the unpolarized 
standards. 

The second method for removing these spectral fringes in¬ 
volves computing the power spectrum of each individual Stokes 
parameter and then filtering unwanted frequencies. We compute 
the power spectrum conventionally as ABS(FFT(q'Mv))^. 

The quv fringe has significant wavelength dependence in the 
800 to lOOOnm region. We choose 817nm, 875nm, 934nm and 
992nm to illustrate the wavelength dependence. We use small 
wavelength intervals of 1000 spectral pixels centered on the 
bandpasses to illustrate the changes with wavelength. Figure [12] 
shows the power spectra for Stokes q, u, and v in all four wave¬ 
length bandpasses. 

There are typically multiple frequency components, partic¬ 
ularly in the Stokes v measurements. For instance, the Stokes 
qu spectra have substantial power at a 1.56 nm ' period in the 


875nm bandpass. The Stokes v spectra are more complex in the 
same bandpass with power at both 1.56 nm * and 0.90 nm“^ pe¬ 
riods. If you model this v power spectrum as a multiplication 
of two cosine functions, you get periods of 0.33 nm ' and 1.23 
nm ' which arise to the sum-frequency of 1.56 nm“^ and the 
difference-frequency of 0.90 nm '. This simple model is shown 
as the solid black curves in every panel of Figure [T2| Since the 
Stokes V data has two retarders mounted in the beam, complex 
interactions are expected. 

The wavelength dependence of the fringe is most readily 
seen as the change in period of the oscillations. The relative 
strength of the fringe at various periods also changes. At 817nm, 
the fringe for Stokes v is largely contained at periods of 1 nm“^ 
and 2 nm ' but with two separate frequency peaks seen at each 
period. By 875nm, the Stokes v power spectra frequencies have 
shifted to shorter frequencies and there are two separate but iden¬ 
tifiable peaks within each main frequency band. At longer wave¬ 
lengths of 934nm, there are two very clear independent frequen¬ 
cies smaller than 1 nm“^ and two much smaller peaks around 1.5 
nm '. The fringe variation in q and u in the 934nm bandpass is 
barely detectable. However, in the 992nm bandpass, the q and u 
variability are clearly detectable again. 

We compared the Fourier filter method with a direct subtrac¬ 
tion of a stellar calibration observation and find good agreement 
where we have sufficient SNR in the Ceres spectrum. However, 
the method if directly subtracting standard star observations un¬ 
der the same optical configuration seems to be robust and to de¬ 
liver higher SNRs when calibration standards are observed at 
very high SNRs. 

Figure [T3] shows the EV Lac quv spectra after subtraction 
of the spectral fringe using two different sets of calibrators. 
We tested the fringe subtraction methodology by using different 
groups of observations to create an average spectrum at much 
higher SNRs. In one group, we simply averaged all observations 
from Table [1] after applying the various flexure and wavelength 
drift corrections to create an average fringe spectrum. For a sec¬ 
ond fringe spectrum, we averaged only the non-brown-dwarf tar¬ 
gets at high SNR observations from Table [1] which excludes the 
LSRJ and 2MASS stars. These two spectra allow us to verify 
consistent results in effectively subtracting this spectral fringe. 

Figure [T3] shows that the resulting fringe subtracted EV Lac 
quv spectra are essentially identical. The SNR is estimated us¬ 
ing the standard deviation in each quv spectrum on both the 
blue and red continuum. The blue continuum runs from 812.2nm 
to 816.3nm while the red continuum runs from 820.9nm to 
824.9nm with each band covering 70 spectral pixels outlined by 
the dashed blue vertical lines in Eigure [13] This method gives 
SNRs from 950 up to 1640 for the red curve and 1310 up to 
1800 for the black curve. As an independent test of the fringe 
subtraction, we apply a high-pass filter to the data by subtracting 
a 9-pixel boxcar smoothed spectrum from the fringe subtracted 
data sets. This 9 pixel smoothing width represents the full width 
10 % max of the delivered optical instrument profile measured 
from arc line exposures. These residual variations after filtering 
are dominated by shot noise and any residual cosmic ray dam¬ 
age. The SNRs of these high pass filtered data sets are in the 
range of 1700 to 2200 showing effective removal of the fringes 
down to typical statistical limits of 0.05% in any individual quv 
spectrum. 

6.1. Slit stepping for High SNR 

An additional method for increasing the SNR is to step the tar¬ 
get along the slit length during an exposure. An example of this 
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Fig. 12. The EV Lac power spectra for the quv fringe in select spectral bandpasses. The central wavelengths are 817nm, 875nm, 934nm and 
992nm with 1000 spectral pixels per computation. Blue shows Stokes q. Green shows Stokes u. Red shows Stokes v. The Stokes v power spectra 
in red show multiple peaks that have significant frequency width. The dashed lines show a simple model with two cosine functions multiplied 
together to give a measure of the single-frequency peak width. The cosine functions have periods of 0.33 nm“' and 1.23 nm“' which arise to the 
sum-frequency of 1.56 nm“* and the difference-frequency of 0.90 nm“* in the power spectrum. 


procedure is shown in Figure [14] The guider offset script was 
set to move the target along the slit in 1 arc second steps dur¬ 
ing an exposure. A mode like this allows the user to achieve 
near-saturation brightness levels on the detector over an order of 
magnitude more spatial pixels without requiring readout of the 
detector. For LRISp on bright targets, the CCD readout time is 
a major limitation. Readout can take up to a minute while the 
integration time to saturation is a small fraction of this time. By 
implementing this slit-stepping mode, the duty-cycle and effi¬ 
ciency of bright target observing campaigns can be kept high 
with integration times substantially longer than the readout time. 

This slit-stepping mode resulted in substantial change in the 
telescope beam footprint as the light passed through the re¬ 
tarders. This changes the spectral fringe. A simple Fourier filter 
as shown above is sufficient to remove spectral fringes. An ex¬ 
ample of a high SNR Stokes q and intensity spectrum for a bright 
held star HD345495 is shown in Figure[T5| After fringe removal, 
the spectra were averaged spectrally by 4 pixels to show SNRs 
above 3000 at 1-point per resolution element spectral sampling. 

We acquired 5 separate spectra of this held star using this 
slit-stepping mode. The guider performance did result in some 
small variation in between steps, but the SNR’s for each individ¬ 
ual polarization measurement were between 1200 and 1500 at 
full spectral sampling. After combining and spectrally averaging 
by a factor of 4, we achieve a hnal SNR of 4500 for a polari- 


metric sensitivity of 0.022% as seen in the blue curve of Figure 
(Bj The higher SNRs achieved by spectral and temporal averag¬ 
ing demonstrates that the SNRs are statistically limited and not 
dominated by some systematic errors. 

6.2. Internal Calibrations 

The coordinate reference frame for linear polarization can 
be verihed with optics m ounted in the internal hlter wheel 
(iGoodrich & CohCTl2003h . The wheel includes two calibration 
polarizers in addition to the QWP used for circular polariza¬ 
tion. The infrared pol arizer is said to have a wa velength range 
of 750nm to 1050nm (IGoodrich & Coh^l2003h . We used this 
calibration optic to perform an independent assessment of the 
cross-talk between linear polarization states introduced by hav¬ 
ing chromatic fast axis orientation variations in the HWR 

We follow a standard method to identify the HWP fast axis 
orientation as well as the overall degree of polarization delivered 
by the calibration optics. The flat field lamps are used to illumi¬ 
nated the fixed calibration IR polarizer. We then take a series of 
exposures while rotating the HWP. In our case, we rotated the 
HWP from -7° to 89° in steps of 2°. This allows us to have high 
angular sampling as well as recording both the minimum and 
maximum intensity transmitted through both beams of the po¬ 
larizing beam splitter. With 49 independent measurements at a 
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Fig. 13. A comparison of the EV Lac quv spectra after fringe sub¬ 
traction using different combinations of stellar and standard observa¬ 
tions. The red curve shows the EV Lac spectra after subtraction of 
all flexure-corrected stellar observations (all stars and calibrators). The 
black curves show the EV Lac spectra when using only bright star and 
calibrator observations (V 2054 Oph, Ceres and unpolarized stars HD 
20630 and HD 174160, no brown dwarfs). The dashed blue lines define 
two regions used for computing SNR statistics as well as continuum 
normalization for the intensity profile. The blue continuum runs from 
812.2nm to 816.3nm while the red continuum runs from 820.9nm to 
824.9nm covering 70 spectral pixels in each band. 
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Fig. 14. An example of the slit-stepping script for achieving high SNR 
measurements. The star was moved in 1 arc second increments along 
the length of the slit 10 times for a total of 11 samples of the local 
seeing and telescope jitter. Each color represents a different exposure. 
One of the two polarized beams had substantially higher throughput. 
Variations along the spatial pixels represents the changing flux on 1 
second timescales. 
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Fig. 15. An example of the q and intensity spectra for a field star using 
the new slit-stepping script. This script allowed us to achieve high SNR 
measurements (>3000). The star was moved in 1 arc second increments 
along the length of the slit following Figure[T4lbefore readout, allowing 
for high duty-cycle observations of bright targets without saturation. 
The blue curve shows an average of 5 repeated exposures on this target 
with SNRs ranging from 1200 to 1500 each. 


range of HWP orientations, we can model the transmitted inten¬ 
sity as a simple function and fit for the calibration parameters. 

At every wavelength, we model the system as a combination 
of an imperfect polarizer and a HWP that rotates the polariza¬ 
tion by some angle. We also include a fit for the period of the 
modulation function with HWP angle as a test on the accuracy 
of the model. Imperfections in the assumptions about the polar¬ 
izer behavior as well as errors in rotation stage encoders values 
can manifest as deviations. The functional form of the intensity 
modulation with HWP angle is represented as an unmodulated 
intensity (/q) constant plus a modulated intensity (/,„) amplitude 
times the modulation function itself. For an ideal HWP with a 
varying fast axis orientation, the modulation is simply a Cosine 
function with an unknown orientation (6). The angular period of 
the polarization modulation (P) should be 90° for a perfect HWP. 
The functional form is thus: 


/ = /o + /„CC5(27r^-^) (9) 

As a test of the method, we perform fits both with a fixed 
P - 90° and a variable period. The results of the functional fit at 
every wavelength are shown in Figure[T6] We included fits where 
P is both fixed and allowed to vary. The fast axis orientation 
variation is seen in color and varies by about 2.5° from 780nm 
to 950nm. In the 950nm to 1050nm the fits begin to diverge. 
The angular modulation period deviates from 90° and the two 
orthogonally polarized beams do not give the same location for 
the HWP fast axis. The assumption of the calibration polarizer 
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ternal calibration optics and can be used without relying on stan¬ 
dard stars. 

During our observing run, we obtained permission to open 
the Keck dome before sunset. We observed the Zenith with 
LRISp before and during twilight on both August 22nd and 23rd. 
During this time, the sky degree of polarization changes quickly 
and takes on spectral features in atmospheric absorption bands. 
The shadow from the earth propagates upward through the atmo¬ 
spheric column during twilight. On August 22nd we observed 
the Zenith sky polarization with the telescope pointed at an az¬ 
imuth of 270° and then 180°. On August 23rd, we observed the 
Zenith at azimuths of 360° and then 090°. Figure [TTI shows the 
detected brightness for two complete polarimetric data sets on 
August 22nd. 


Fig. 16. The degree of polarization for the IR calibration polarizer and 
the half-wave plate fast axis orientation. The measured degree of po¬ 
larization (DoP) for the polarizer is shown on the right and y-axis. The 
DoP is above 95% for wavelengths shorter than 900nm. At longer wave¬ 
lengths the DoP falls to below 20% at lOOOnm. The fast axis orienta¬ 
tion for the HWP is shown as the colored curves with the left-hand y 
axis. The colors correspond to the two orthogonally polarized beams 
recorded on the detector produced by the polarizing beam splitter. Each 
colored fast axis location has two separate curves corresponding to fit¬ 
ting a function with either a fixed angular modulation period {P = 90°) 
or with a variable period (P). The fast axis orientation changes by about 
2.5° from 780nm to 950nm. Beyond 950nm, the calibration polarizer 
has low DoP and the fast axis orientations measured in the two beams 
starts to diverge. Additionally, the computed modulation period diverges 
from the nominal angle of P = 90°. Likely the low DoP and other po¬ 
larization artifacts from the polarizer itself cause the measurement tech¬ 
nique to give errors. See text for details. 


being simply represented as a fractional polarization at a fixed 
angle and the behavior of the HWP appears to break down. 

However, this testing does show that the expected rotation of 
the linear polarization reference frame caused by deviations in 
HWP properties with wavelength is well controlled. The rotation 
of the plane of polarization for a science target should be easily 
compensated by use of a linear pola rization standard star as is 
common for calibrat ions of LRISp (iGoodrich & Coh^l2003l : 
iGoodrich et al.llToO^ . 


7. Daytime Sky Calibration Testing: Cross-talk 

As an independent test of the instrument calibration, we use 
the daytime sky polarization. The quarter wave plate used to 
measure circular polarization is mounted in the calibration filter 
wheel. As such, there are no calibration optics mounted in front 
of the instrument in the 2-retarder configuration. This presents a 
calibration challenge. Recently, we have been developing meth¬ 
ods to use daytime sky polarization as a bright, highly polar¬ 
ized calibration source with a well kno wn angle of polarization 
(AGP ) to derive telescope properties (iHarrington et al.ll201lL 
I2OIOI) . We have developed algorithms to use this daytime sky po¬ 
larization to compute the cross-talk introduced in the instrument 
while illuminating the entire optical train much more similar to 
star light. This method was successful in calibrating the AEOS 
telescope and the HiVIS spectropolarimeter where the linear to 
circular cross talk was 100% at some wavelengths and telescope 
pointings. This method provides an independent check of the in- 



Fig. 17. Daytime and twilight sky intensity measured with LRISp as 
the sun set on August 22nd. The first data set was recorded with the 
telescope at an azimuth of 270 shown in solid lines. The second data set 
was recorded at an azimuth of 180 shown in dashed lines. Each color 
shows a different modulation state. Note there is at least an order of 
magnitude change in brightness between spectra just due to the high 
degree of polarization in the daytime sky at the Zenith. The brightness 
dropped by 3 orders of magnitude between the start and end of the test 
as the sun set. 


At the Zenith on a mountain site with the sun near the hori¬ 
zon, we expect to measure degree of polarizations above 60% 
as measured by all-sk y imaging polarimeters common in the 


atmospheric sciences (ISwindlel 20141 ISwindle & KuhnI 

2014 

Dahlberg et al.l2009l.l201 lUShaw et alJ2010l:lPust & Shaw 

2006 

2007. 2 OO 9 I). Figurell8lshows the Stokes quv spectra and the as- 


sociated Degree of Polarization (DoP) measured with LRISp on 
August 22nd. 


We measured 65% to 70% DoP in the first data set and 80% 
in the second set. This slight increase in DoP as the sun reaches 
the horizon is expected from standard Rayeligh scattering theory 
on a clear day. 

We can see that the Stokes v measurements are small but 
nonzero. This is expected as there is a known misalignment be¬ 
tween the retarders and the polarizing beam splitter orientation 
in addition to th e expected chromatic change in the QWP fast 
axis orientation (IGoodrich et alJ[T99-5h . Note that the QWP is 
fixed in the calibration wheel mount and optimization must be 
done manually in the present instrument configuration. The day¬ 
time sky has no circular polarization measured to limits better 
than 1% (ISwindlell2014t [dwindle & Kuhnll2014l) . Our measure¬ 
ments show that the linear to circular cross-talk is roughly a few 
percent of the incoming linear polarization signal. 
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Fig. 18. The Stokes quv spectra and the corresponding degree of po¬ 
larization (DoP) measured with LRISp on August 22nd. The solid lines 
show the first data set with the sun above the horizon and the telescope 
at an azimuth of 270. The dashed lines show the second data set with 
the telescope at an azimuth of 180 as the sun was setting. The Stokes 
V measurement (dashed red line) has the lowest flux levels and the cor¬ 
responding highest noise at longer wavelengths. Clear linear-to-circular 
cross talk is seen as the nonzero Stokes v component of roughly 5% 
changing sign as the instrument rotates. 


We also note that the angular relationship was verified to be 
preser ved following our algorithms outlined in iHarrington et al.l 
(1201 ih . On both August 22nd and 23rd we observed the Zenith 
with the telescope at cardinal pointings of North, East, South 
and West. The angle between the measured Stokes vector and 
the theoretical Rayleigh sky vector was preserved to better than 
5° as the instrument rotated on both days. 

The high measured degree of polarization for the day¬ 
time sky shown above at long wavelengths also suggests effi¬ 
cient modulation and polarizing beamsplitter functionality even 
though the internal calibration polarizer was unable to accurately 
verify long wavelength performance from the degraded internal 
polarizer performance. 

7.1. Intensity to Polarization Cross-talk 

Another critical test for any spectropolarimeter is to measure the 
level of cross-talk between the intensity spectrum and the po¬ 
larization spectra. There are two separate types of intensity to 
polarization artifacts. The instrumental induced polarization is 
the first type. This real polarization signal is generally caused by 
the optics. This is often called the zero point calibration and is 
typically measured with an unpolarized standard stars. This po¬ 
larization should not reflect the incoming stellar line spectrum 
and is usually a smooth continuum function. As with most spec¬ 
trographs our unpolarized standard measurements show ~1% in¬ 
strumental polarization as relatively smooth functions of wave¬ 
length. 

A problem with defocusing the telescope is that the defo- 
cused beam as masked by the slit creates unstable zero point 
instrument polarization when combined with the poor guiding 
performance of the Keck slit guider. Defocusing as a means for 
increasing the dynamic range has several advantages for line 
polarization studies where the continuum is simply fit and sub¬ 
tracted. However, the relatively unstable masking of the incom¬ 
ing beam breaks the circular symmetry of optical path and re¬ 


duces the ability for the instrument to measure continuum polar¬ 
ization. Care must be taken when defocusing the telescope as the 
continuum stability is reduced. In addition, the slit guider per¬ 
formance degrades at high airmass and continuum polarization 
is less stable when tracking targets at high elevations. 

A second kind of artifact is caused by instrument and data re¬ 
duction artifacts. This kind of intensity to polarization cross talk 
involves having the incoming stellar spectrum spuriously present 
in the polarization spectrum. Imperfections in detector linearity, 
background subtraction and other effects begin to cause limita¬ 
tions at high sensitivity levels c.f. lKelle^ (Il996l) . In unpolarized 
standard stars, in addition to the instrumental continuum polar¬ 
ization, the stellar line spectrum along with it’s first and second 
wavelength derivatives can be spuriously present. Careful cor¬ 
rection for scattered light, ghost images, wavelength drifts, etc 
will only reduce this kind of I to QUV cross talk below some 
detection threshold. The dual-beam configuration in addition to 
beam-swapping during data reduction does minimize several ar¬ 
tifacts to first order. However, several second order effects are 
not completely removed and can easily dominate the quv spectra 
at high SNR levels. 
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-50 0 50 
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Fig. 19. The cross correlation functions computed with the IDL 
C_CORRELATE routine. The black lines show the correlation func¬ 
tions between intensity and the qu spectra for the 5 slit-stepped spectra 
of HD 354495 described above. No substantial peaks are seen above 
the typical fluctuation levels of 0.1. The blue lines show the correlation 
functions when 0.02% of the continuum normalized intensity spectrum 
was added to the individual polarization spectra. These blue correlation 
curves show that a substantial correlation of 0.2 to 0.3 is present at zero 
lag giving rise to our upper limits for intensity to polarization cross-talk. 


For LRISp, we took the five high SNR measurements in slit- 
stepped mode used to create Figure [15] By running intensity- 
to-polarization correlations across many spectral lines, a clear 
signature of intensity to quv spectra can be assessed. Each of the 
five spectra had a SNR of 1200 to 1500 in each q and u Stokes pa¬ 
rameter. A correlation analysis of 1000 spectral pixels between 
856nm and 896nm yielded no significant I to qu cross-talk at lev¬ 
els above the qu SNR. The same test was run on the continuum 
normalized polarization profiles Q — I * q (sometimes differen¬ 
tiated as qc as opposed to q). A test of the routines showed that 
0.2% of the intensity spectrum added in to either qu or QU spec¬ 
trum was easily detectable with our correlation analysis. Figure 
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[T^shows these normalized correlation curves computed with the 
IDL C_CORRELATE function. A correlation of 1.0 would de¬ 
scribe perfect correlation. The black curves show the correlation 
between the continuum normalized intensity spectra and the qu 
spectra. The correlation curves fluctuate around zero with low 
amplitudes of 0.1 or less showing no substantial correlation. The 
blue curves of Eigurejl^show the correlation functions with the 
intensity added to the qu spectra with a 0.2% multiplier. Clear 
peaks are seen with correlations above 0.2 when we simulate in¬ 
tensity to polarization cross-talk at 2x the SNR levels (0.2%). 
This test sets our upper limits on the I to qu cross-talk for our 
instrument configuration to <0.2%. 

8. Comparison: LRISp Blue Channel 

The blue channel of the polarimeter has some substantial differ¬ 
ences in sensitivity to polarimetric errors when using common 
spectral extraction algorithms. We used the 680 dichroic and the 
300/5000 grism during this campaign. The full-width half-max 
of arc line fits is between 6.0 and 7.5 spectral pixels showing 
substantial oversampling. There are slight differences of roughly 
0.1 to 0.3 pixels EWHM between the two orthogonally polarized 
beams, showing a small but detectable difference in optical reso¬ 
lution between the orthogonally polarized beams. The delivered 
spectral resolution is R~450 with 6 points per resolution element 
at 360nm. The resolution rises to R~790 at 7 point per EWHM 
sampling at 775nm. 



Fig. 20. A comparison between the pixel coordinate grid tilt geomet¬ 
rical calibrations derived using the arc lamp fits. Black symbols show 
the computed wavelength center of several arc lines at selected spatial 
pixels on the CCD. The red symbols show the median spatial location 
of all arc line wavelength centers for the two orthogonally polarized 
beams of the red channel. Because the average wavelength solution in 
fixed ccd pixel coordinates was used, tilts between the orthogonally po¬ 
larized beams are seen between a fixed pixel grid and the tilted arc lines. 
This coordinate geometry must be corrected in the red channel by asso¬ 
ciating each spatial pixel with the correct wavelength derived from the 
arc lines as described above. The blue channel has a greatly reduced tilt, 
as shown by the blue symbols for the median arc line location across the 
blue channel CCD. 

Unlike the red channel, the tilt of the monochromatic slit im¬ 
age is less than 0.1 pixels across 240 spatial pixels extracted. 
Tilted coordinate geometries can be overcome using up sam¬ 
pling and instrument calibrations, or even more complex instru¬ 
ment profile deconvolutions. However, with such small spectral 
tilt in this blue channel, the wavelength errors and polarimetric 


artifacts are greatly reduced even when using simple extraction 
algorithms. The red and blue channel spectral tilts are shown 
in Eigure The red tilt is over 2 pixels across the slit image 
while the blue channel is more than an order of magnitude less. 
Arc line exposures are used to determine how each wavelength 
falls across the CCD pixel grid. The arc line central wavelength 
is mapped at every spatial and spectral pixel for every arc lamp 
line, shown as the black symbols in Eigure |20] only for the red 
channel. The median spatial locations of all arc lamp lines along 
the CCD pixel grid are shown in red for the LRISp red channel. 
The same median location for the blue channel is shown in blue. 

The blue channel ccd is much thinner than the red channels 
deep depletion pixels. The corresponding cosmic ray hit rate is 
greatly reduced. The optimal extraction algorithms outlined for 
the red channel are always effective for rejecting cosmic ray hits 
as well as filtering some detector cosmetic issues. However, there 
is less need for this algorithm for modest exposure times. 
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Fig. 21. The spatial profiles computed for EV Lac. The average spatial 
profile over all spectral pixels is shown for all 6 polarimetric exposures 
for the red channel (plotted in red) and the blue channel (plotted in blue). 
The blue channel has a ghost image and substantially wider scattered 
light in the wings of the spatial profile. 


Careful scattered light and background subtraction is more 
important in the blue channel. There is a ghost image and sub¬ 
stantially more scattered light width to the blue channel spa¬ 
tial profile. A comparison between the blue and red channels is 
shown in Figure |2T] The blue channel takes roughly double the 
number of spatial pixels to drop below the 1% core flux level. 
Typical background subtraction algorithms will include some 
fraction of the stellar flux depending on the reduction choices 
made. This is easily subtracted along with night sky backgrounds 
provided this scattered light profile is known and properly ac¬ 
counted in the reduction parameter choices. Night sky lines must 
be assumed to sit on top of a scattered light profile that extends 
substantially away from the core region, and scattered light back¬ 
grounds must be measured from an appropriate spatial distance 
away from the stellar core. 

The curve of growth is used to determine the optimal spatial 
extraction width when considering different noise source. We de¬ 
fine the curve of growth as the computation of how much flux is 
included in the spectrum when setting progressively wider spa¬ 
tial profiles divided by the total detected flux. For instance, if the 
user picks a blue channel half width of 10 for a total of 21 spatial 
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pixels in the blue extracted spectrum, the user is only including 
90% of the total flux. If the user would select a half width of 39 
spatial pixels, 99.9% of the light would be included. In the red 
channel, care must be taken as cosmic ray damaged pixels are in¬ 
completely corrected and this noise source can dominate errors 
from incomplete flux inclusion. Thankfully the scattered light is 
reduced in the red channel and correspondingly less spatial pix¬ 
els are required to gather the majority of the detected stellar flux. 

Spectral fringe is not seen in the blue channel at these low 
spectral resolutions. An example high sensitivity polarimetric 
data set for EV Lac is shown in Figure |22] The signal to noise 
levels achieved are above 2000 per spectral pixel at full spectral 
sampling as estimated from the quv spectra. The quv spectra are 
dominated by photon statistics across the entire 400nm to 700nm 
range. The 680nm dichroic cutoff used did have some leakage 
beyond the cutoff wavelength, and spectra were extracted out to 
775nm. In the EV Lac spectra, a small spectral region was satu¬ 
rated and was subsequently not plotted. 



Fig. 22. The quv and intensity spectra for EV Lac recorded on August 
22nd 2012. The signal to noise ratios are estimated to peak above 2000 
as seen in the quv noise levels. 


The blue channel showed similar cross-talk levels to the 
red channel given the standard retarder alignment and operat¬ 
ing angles. The daytime sky polarization observations were per¬ 
formed both on red and blue channels. The measured sky degree 
of polarization is between 50% and 80% for the entire 360nm 
to 770nm range. Figure |23] shows the daytime sky polarization 
measurements. Cross talk between linear and circular states is 
seen at similar levels to the red channel. 

9. Summary 

We have developed a data reduction pipeline that can calibrate 
the LRISp spectropolarimeter and achieve signal-to-noise ra¬ 
tios above 2000 limited by photon statistics at full resolution 
and sampling. Spectral binning and temporal averaging achieved 
SNRs of 4500 without obvious visible instrumental artifacts. 



Fig. 23. The daytime sky polarization measured roughly 30 minutes be¬ 
fore sundown on August 22nd 2012. Blue shows q, yellow shows u and 
red shows v. The total degree of polarization (DoP) is shown in black. 
The 680nm dichroic cutoff greatly reduces the detected sky brightness 
for longer wavelengths. Poor background subtraction and the blue chan¬ 
nel scattered light profiles contribute to the change in polarization be¬ 
havior given standard reduction algorithms. 


We have demonstrated algorithms for overcoming several in¬ 
strument artifacts and tested them on both red and blue chan¬ 
nels. The major polarimetric limitations proved to be wavelength 
drifts from flexure, spectral fringes from the retarders and ef¬ 
fective removal of cosmic ray contamination in the red channel 
deep depletion CCDs. With these calibrations, we have success¬ 
fully reproduced magnetic held measurements in atomic lines 
of M dwarfs such as EV Lac and found new magnetic signa¬ 
tures in molecular bands of M dwarfs and brown dwarfs. Sensi¬ 
tive comparisons of low resolution spectropolarimetric data with 
new magnetic held models such as Figure |24] can now begin. 
With magnetic fields producing signals at the >0.2% levels, we 
must achieve shot-noise limited performance to sensitivity levels 
substantially below this level as demonstrated here. 

The instrument flexure in addition to beam wobble from ro¬ 
tating retarders introduce wavelength instabilities. The drifts are 
a substantial fraction of a pixel within and between modulated 
polarimetric exposures. These systematic errors produce spuri¬ 
ous instrumental artifacts that are proportional to the first and 
second derivative of the spectral intensity with wavelength, even 
if the tilted coordinate geometry is calibrated using standard 
arc lamp exposures. These artifacts exactly resemble signatures 
from stellar magnetic fields and must be effectively suppressed 
to achieve accurate science results. Using correlation techniques 
with telluric lines and atmospheric sky-glow lines we can effec¬ 
tively track and remove these wavelength drifts. 

The spectral fringes introduced by the retarders have ampli¬ 
tudes of 0.2% in Stokes qu and over 0.5% in Stokes v in the red 
channel. These fringes display wavelength dependent behavior 
in both amplitude and frequency for individual Stokes param¬ 
eters. Simple Fourier Alters can remove the fringes provided a 
bandpass of 50nm or less is used. The Fourier Altering method 
gives very similar results to subtracting a calibration standard 
star observation. However, calibration of fringes by subtraction 
of standard star observations does give better results. 

Cosmic ray hits are present in more than half the spectral pix¬ 
els for each Stokes parameter in a typical long exposures of the 
red channel on a faint target. Optimal spectral extraction tech- 
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niques are an effective filter of this noise source. An estimate 
of the local spatial prohle computed using adjacent wavelengths 
gives information about the expected prohle at any individual 
wavelength. This local spatial prohle is used in an iterative loop 
to identify and reject cosmic ray hits and to correct the damaged 
pixels using the best-estimate of the spatial prohle shifted and 
scaled to each individual extracted wavelength. 
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Fig. 24. The model Stokes quv signatures induced by a stellar magnetic 
held. Our brown dwarf target LSRJ was modeled under a range of mag¬ 
netic held orientations, strengths and hlling factors to show examples 
of the quv morphology variation with magnetic held properties. Typical 
kilogauss helds as inferred from radio and optical measurements pro¬ 
duce detectable signatures in the CrH band here with amplitudes of less 
than 0.4%. 


We have applied our daytime sky based polarization calibra¬ 
tion techniques to LRISp to derive estimates of the instrument 
quv cross-talk and polarization reference frame. There is a few 
percent linear to circular cross talk present as seen by observed 
Stokes V when looking at a ~100% linearly polarized daytime 
sky in both blue and red channels. This is expected given the 
wavelength dependent properties of the retarders and small mis¬ 
alignments in the retarder rotation stages. We also assessed the 
intensity to qu cross-talk on high SNR observations and derive 
upper limits. The red channel I to qu and QU cross-talk is below 
0 .2% and limited by shot noise. 

This pipeline will allow high SNR use of LRISp and Keck 
on faint targets now that the major sources of instrumental error 
have been identihed and suppressed. 
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